Spin injection and spin accumulation in all-metal mesoscopic spin valves 
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We study the electrical injection and detection of spin accumulation in lateral ferromagnetic metal- 
nonmagnetic metal-ferromagnetic metal (F/N/F) spin valve devices with transparent interfaces. 
Different ferromagnetic metals, permalloy (Py), cobalt (Co) and nickel (Ni), are used as electrical 
spin injectors and detectors. For the nonmagnetic metal both aluminium (Al) and copper (Cu) are 
used. Our multi-terminal geometry allows us to experimentally separate the spin valve effect from 
other magneto resistance signals such as the anomalous magneto resistance (AMR) and Hall effects. 
We find that the AMR contribution of the ferromagnetic contacts can dominate the amplitude of the 
spin valve effect, making it impossible to observe the spin valve effect in a 'conventional' measurement 
geometry. In a 'non local' spin valve measurement we are able to completely isolate the spin valve 
signal and observe clear spin accumulation signals at T = 4.2 K as well as at room temperature 
(RT). For aluminum we obtain spin relaxation lengths (\ s f) of f.2 fim and 600 nm at T = 4.2 K and 
RT respectively, whereas for copper we obtain 1.0 /xm and 350 nm. The spin relaxation times r s f in 
Al and Cu are compared with theory and results obtained from giant magneto resistance (GMR), 
conduction electron spin resonance (CESR), anti-weak localization and superconducting tunneling 
experiments. The spin valve signals generated by the Py electrodes (ayAf = 0.5 [1.2] nm at RT 
[T = 4.2 K]) are larger than the Co electrodes (a F \ F = 0.3 [0.7] nm at RT [T = 4.2 K]), whereas 
for Ni (of^f < 0.3 nm at RT and T = 4.2 K) no spin signal is observed. These values are compared 
to the results obtained from GMR experiments. 

PACS numbers: 72.25.Ba, 72.25. Hg, 72.25.Mk, 72.25.Rb 



I. INTRODUCTION 

Spintronics is a rapidly emerging field in which one 
tries to study or make esaLicit use of the spin degree of 
freedom of the electron.ElcrH Sofar, the most well known 
examples of spintronics are the fiiwt magneto resistance 
(GMR) of metallic multi-layersclBB and tunneling mai 
neto resistance (TMR) of magnetic tunnel junctionsQI 
Injection of hot electrons w 1 eV above the Fermi en- 
ergy Ep in Co/Cu (multi) layers have shown a significant 
spin filtering effect, enabling transistoij-|f-wictionality and 
ballistic electron magnetic microscopy£H13 Recent exper- 
iments have shown the ability of spin polarized currents 
to initiate a (local) magnetization reversal in i.hjji .fnjra- 
magnetic wires and Co/Cu multi-layer pillarsO£3tj'oli3. 
A new direction is emerging, where one actually wants 
to inject spin currents, transfer and manipulate the spin 
information at the Fermi-level, and detect the resulting 
spin polarization in nonmagnetic metals and semiconduc- 
tors. Because of spin-orbit interaction, the electron spin 
can be flipped and consequently a spin polarized current 
will have a finite lifetime. For this reason it is necessary 
to study spin transport in systems, where the 'time of 
flight' of the electrons between the injector and detec- 
tor is shorter than the spin relaxation time. A first and 
successful attempt to electrically inject and detect spins 
in metals dates back to 1985 when Johnson and Silsbee 
successfully demonstrated spin accumulation in a su 
crystal aluminium bar up to temperatures of 77 KE!; 
In their pioneering experiments they were able to observe 
spin precession of the induced non-equilibrium magneti- 
zation, made possible by the long spin relaxation lengths 



A s / > 50 (im. In (diffusive) thin metallic films however, 
the spin relaxation length corresponds to typical length 
scales of 1 /xm. We use a lateral mesosapie spin valve, 
to access and probe this length scaleEJOEIH. We note 
that a similar experiment using planar spin valves has 
been reported in Ref. [2^. 

In Sec. II a review of the basic model for spin trans- 
port in the diffusive transport regime is given, whereas 
in Sec. Ill this model is applied to our multi-terminal 
device geometry. A multi-terminal resistor model of spin 
injection and detection is presented in Sec. IV in order to 
elucidate the principles behind the reduction of the po- 
larization of the spin current at a transparent F/TOnter- 
face, also referred to as "conductivity mismatch" r3 The 
sample fabrication process and measurement geometry 
are described in Sec. [v|. Spin accumulation measure- 
ments in a 'conventional' and 'non-local' geometry for 
Py/Cu/Py and Py /Al/ Py spin valves will be presented 
in Sec. [Vjj and Sec. VII, whereas spin accumulation mea- 



surements on Co/C u/Co and Ni/Cu/Ni spin valves will 
be presente d in S ec. |VIII|. In Sec. IX the obtained results 
of Sees. VI VII and [VIIl| are analyzed using the model 



for spin transport in the diffusive regime and the results 
are compared to GMR, CESR, anti-weak localization and 
superconducting tunneling experiments. 



II. THEORY OF SPIN INJECTION AND 
ACCUMULATION 

In general, electron transport through a diffusive chan- 
nel is a result of a difference in the (electro-)chemical po- 
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tential of two connected electron reservoirs. An electron 
reservoir is an electron bath in full thermal equilibrium. 
The chemical potential fi c h is by definition the energy 
needed to add one electron to the system, usually set 
to zero at the Fermi energy (this convention is adapted 
throughout this text), and accounts for the kinetic en- 
ergy of the electrons. In the linear response regime, i.e. 
for small deviations from equilibrium (\eV\ < kT), the 
chemical potential equals the excess electron density n 
divided by the density of states at the Fermi energy, 
Mc/i = n/N(Ep). In addition an electron may also have 
an potential energy, e.g. due to the presence of an electric 
field E. The additional potential energy for a reservoir at 
potential V should be added to fi c h in order to obtain the 
electrochemical potential (in the absence of a magnetic 
field): 



fi = fi ch - eV , 



(1) 



where e denotes the absolute value of the electron charge. 

From eq. jjj it is clear that a gradient of //, the driv- 
ing force of electron transport, can result from either a 
spatial varying electron density Vn or an electric field 
E = — VV. Since /x fully characterizes the reservoir one 
is free to describe transport either in terms of diffusion 
(E = 0, Vn ^ 0) or in terms of electron drift (E ^ 0, 
Vn = 0). In the drift picture the whole Fermi sea has 
to be taken into account and consequently one has to 
maintain a constant electron density everywhere by im- 
posing: Vn = 0. We use the diffusive picture where only 
the energy range A/i, the difference in the electrochemi- 
cal potential between the two reservoirs, is important to 
describe transport. Both approaches (drift and diffusion) 
are equivalent in the linear regime and are related to each 
other via the Einstein relation: 



a = e 2 N{E F )D 



(2) 



where a is the conductivity and D the diffusion constant. 

We focus on the diffusive transport regime, which ap- 
plies when the mean free path l e is shorter than the device 
dimensions. The description of electrical transport in a 
ferromagnet in terms of a two-current (spin-up and spin- 
down) model dates back to Fert and CampbclO. Van 
Son et alr3 have extended the model to describe trans- 
port through fcrromagnct-nonmagnetic metal interfaces. 
A firm theoretical underpinning, based on the Boltz- 
mann_transport equation has been given by Valet and 
Fert .Ell They have applied the model to describe the ef- 
fects of spin accumulation and spin dependent scattering 
on the GMR effect in magnetic multilayers. This stan- 
dard model allows for a detailed quantitative analysis of 
the experimental results. 

An alternative model, based on thermodynamic con- 
siderations, has been put forward and applied by Johnson 
and Silsbee (JS)E3. In principle both models describe the 
same physics, and should therefore be equivalent. How- 
ever, the JS model has a drawback in that it does not 



allow a direct calculation of the spin polarization of the 
current (77 in Refs. [f6| , p8 29 3(]), whereas in the standard 
model all measurable quantities can be directly, related 
to the parameters of the experimental system.EJ'EjEj 

The transport in a ferromagnet is described by spin 
dependent conductivities: 



cr T = V|g 2 Z3|, with _Dj = —Vp^lel (3) 
a i = N ie 2 D u with D i = ^v Fl l el , (4) 

where iV-j-j denotes the spin dependent density of 
states (DOS) at the Fermi energy (Ep), and £>+ j_ the 
spin dependent diffusion constants, expressed in the spin 
dependent Fermi velocities vft,1j an d electron mean free 
paths let, l- Throughout this paper our notation is f for 
the majority spin direction and J, for the minority spin 
direction. Note that the spin dependence of the conduc- 
tivities is determined by both density of states and diffu- 
sion constants. This should be contrasted with magnetic 
F/I/F or F/I/N tunnel junctions, where the spin polar- 
ization of the tunneling elentmiis, is determined by the 
spin-dependent (local) DOS.QSl-J Also in a typical fer- 
romagnet several bands (which generally have different 
spin dependent density of states and effective masses) 
contribute to the transport. However, provided that the 
elastic scattering time and the interband scattering times 
are shorter than the spin flip times (which is usually the 
case) the transport can still be described in terms of well 
defined spin up and spin down conductivities. 

Because the spin up and spin down conductivities are 
different, the current in the bulk ferromagnet will be dis- 
tributed accordingly over the two spin channels: 



J] 

J I 



e dx 
e dx 



(5) 
(6) 



where are the spin up and spin down current den- 
sities. According to Eqs. [| and || the current flowing in 
a bulk ferromagnet is spin polarized, with a polarization 
given by: 



Qf = 



(7) 



The next step is the introduction of spin flip processes, 
described by a spin flip time r-fj, for the average time 
to flip an up-spin to a down-spin, and tj,j for the re- 
verse process. The detailed balance principle imposes 
that N^/t^i = -V|/t||, so that in equilibrium no net 
spin scattering takes place. As pointed out already, usu- 
ally these spin flip times are larger than the momentum 
scattering time r e = l e /vp. The transport can then be 
described in terms of the parallel diffusion of the two 
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spin species, where the densities are controlled by spin 
flip processes. It should be noted howeitet-that in partic- 
ular in ferromagnets (e.g. permalloyE3'E3E3) the spin flip 
times may become comparable to the momentum scat- 
tering time. In-.this.case an (additional) spin-mixing re- 
sistance arisesoHE3, which we will not discuss further 
here. 

The effect of the spin flip processes can now be de- 
scribed by the following equation (assuming diffusion in 
one dimension only): 



D 



dx 2 



(8) 




+ L/2 < 
IV + t 



■ L/2 



VI 



where D = D^D l (N^+N l )/(N^D- l +N l D L ) is the spin 
averaged diffusion constant, and the spin relaxation time 
r s f is given by: l/r s f = 1/Vf j + 1/tjj. We note that T s f 
represents the timescale over which the non-equilibrium 
spin accumulation (/if — /ij.) decays and therefore is equal 
to the spin lattice relaxation time T\ used in the Bloch 
equations: T s f = T\ .cj Using the requirement of current 
conservation, the general solution of eq. || for a uniform 
fcrromagnet or nonmagnetic wire is now given by: 



C D 
/if = A + Bx^ exp(—x/X s f) H exp(x/\ s f) (9) 

<7f <7f 

C D 
p\ = A + Bx exp(—x/X s f) exp(x/\ 8 f)().0) 



where we have introduced the spin relaxation length 
A s / = y/ Dr s f ■ The coefficients A,B,C, and D are deter- 
mined by the boundary conditions imposed at the junc- 
tions where the wires are coupled to other wires. In the 
absence of an interface resistance and spin flip scattering 
at the interfaces, the boundary conditions are: 1) conti- 
nuity of /if, pi at the interface, and 2) conservation of 
spin-up and spin-down currents fa, j± across the inter- 
face. 



III. SPIN ACCUMULATION IN 
MULTI-TERMINAL SPIN VALVE STRUCTURES 



We will now apply the model of spin injection to a 
non local geometry, which reflects our measurement and 
device geometry, see Fig.jlja and Fig||c. 

In our (I-dimensional) geometry we can identify 6 dif- 
ferent regions for which Eqs. ^ and [lO] have to be solved 
according to their boundary conditions at the interface. 
The geometry is schematically shown in Fig^TJo, where 
the 6 different regions are marked with roman letters I 
to VI. According to Eq. |9| the equations for the spin- 
up electrochemical potentials in these regions, assuming 
parallel magnetization of the ferromagnetic regions, read: 



FIG. 1: (a) Schematic representation of the multi-terminal 
spin valve device. Regions 7 and VI denote the injecting (Fx) 
and detecting (F2) ferromagnetic contacts, whereas regions 
77 to V denote the four arms of a normal metal cross (TV) 
placed in between the two ferromagnets. A spin polarized 
current is injected from region 7 into region 77 and extracted 
at region IV. (b) Diagram of the electrochemical potential 
solutions (Eqs. |^ and |l^) in each of the six regions of the 
multi-terminal spin valve. The nodes represent the origins of 
the coordinate axis in the 6 regions, the arrows indicate the 
(chosen) direction of the positive x-coordinate. Regions 77 
and 777 have a finite length of half the Py electrode spacing 
L. The other regions are semi-infinite. 
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where we have written <7f = <jf(1 + ap)j2 and 
A, B,C, D, E, F,G, H and K are 9 unknown constants. 
The equations for the spin-down electrochemical poten- 
tial in the six regions of fig. can be found by putting a 
minus sign in front of the constants C, D, E, F, H, K, G 
and ap in Eqs. I to VI. Constant B is the most valuable 
to extract from this set of equations, for it gives directly 
the difference between the electrochemical potential mea- 
sured with a normal metal probe at the center of the 
nonmagnetic metal cross in fig.[j]a and the electrochemical 
potential measured with a ferromagnetic voltage probe at 
the F/N interface of region V and VI. For X s f >> L i.e. 
no spin relaxation in the nonmagnetic metal of regions II 
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and V, the ferromagnetic voltage probe effectively probes 
the electrochemical potential difference between spin-up 
and spin-down electrons at center of the nonmagnetic 
metal cross. Solving the Eqs. I to VI by taking the 
continuity of the spin-up and spin-down electrochemi- 
cal potentials and the conservation of spin-up and spin 
down-currents at the 3 nodes of Fig. [l]b, one obtains: 



B 



2 \N_ e -L/2\ N 
e ctjv 



2(M + l)[Msinh(L/2\ N ) + cosh(L/2X N )} ' 

(11) 



where M = (<jfXn /0n^f)(1 — and L is the length 
of the nonmagnetic metal strip in between the ferromag- 
netic electrodes. The magnitude of the spin accumula- 
tion at the F/N interface of region V and VI is given by: 
MT - Ml = B/a F . 

In the situation where the ferromagnets have an anti- 
parallel magnetization alignment, the constant B of Eq. 
|ll| gets a minus sign in front . Upon changing from par- 
allel to anti-parallel magnetization configuration (a spin 
valve measurement) a difference of Afi = 2B will be de- 
tected in electrochemical potential between the normal 
metal and ferromagnetic voltage probe. This leads to the 
definition of the so-called spin-coupled or spin-dependent 



resistance AR 



2B 
-ejS 



, where S is the cross-sectional 



area of the nonmagnetic strip: 



AR 



a 2 XM^ e -L/2X N 

r crjvo 



(M + l)[Msinh(L/2\ N ) + cosh(L/2X N )} 



(12) 



Equation |l2j shows that for Ajy << L, the magnitude of 
the spin signal AR will decay exponentially as a function 
of L. In the opposite limit, Xf << L « Xn the spin 
signal AR has a 1/L dependence. In this limit and under 
the constraint that ML/2 X n >> 1, we can write Eq. |l2| 
as: 



AR = 



2a 2 F X 2 N 



M[M + l)a N SL 



(13) 



In the situation where there are no spin flip events in 
the normal metal (Xn — oc) we find that we can write 
eq. Il3l in an even more simple form: 



AR = 



2a 2 F X F /a 



(1 - a F ) 2 SL/o N 



(14) 



The important point to notice is that Eq. [l4| clearly 
shows that even in the situation when there are no spin 
flip processes in the normal metal, the spin signal AR is 
reduced with increasing L. The reason is that the spin de- 
pendent resistance (Xf/ct f S) of the injecting and detect- 
ing ferromagnets remains constant for the two spin chan- 
nels, whereas the spin independent resistance (L/ctnS) of 



the nonmagnetic metal in between the two ferromagnets 
increases linearly with L. In both nonmagnetic metal 
regions II and V (Fig. the spin currents have to tra- 
verse a total resistance path over a length X F + L/2 and 
therefore the polarization of the current flowing through 
these regions will decrease linearly with L and hence the 
spin signal AR. Note that in the regions V and VI no 
net current is flowing as the opposite flowing spin-up and 
spin-down currents are equal in magnitude. 

Using Eqs. |[ ^ and / we can calculate the current 
polarization at the interface of the current injecting con- 



tact, defined as P 



We obtain: 



P : 



CtF 



Me L/2\ N + 2cosh(L/2X 



N) 



2(M + l)[Msinh(L/2X N ) + cosh(L/2X N )} 



(15) 



In the limit that L » Xn we obtain_the polarization 
of the current at a single F/N interface:^ 



aj? 
M + 1 



(16) 



Again, Eq. |l6| shows a reduction of the polarization 
of the current at the F/N interface, when the spin de- 
pendent resistance (Xf / <JfS) is much smaller that the 
spin independent resistance (Xn/&nS) of the nonmag- 
netic metal. This situation becomes progressively worse 
for a semiconductor as ctat is reduced by a factor of 100 
or more and.. has become known as the "conductivity 
mismatch" .BE 

Finally we note that the spin signal /\R Conv can a l so 
be calculated for a conventional measurement geometry, 
see Fig. [|b, writing down similar equations and bound- 
ary conditions as we have done for the non local geometry 
(Eqs. / to VI). We find: 



AR 



Conv 



= 2AR 



(17) 



Equation [l?] shows that the magnitude of the spin 
valve signal measured with a conventional geometry is 
increased with a factor two as compared to the non local 
spin valve geometry (see also Ref. Bfl Eq. 45). 



IV. RESISTOR MODEL OF MULTI-TERMINAL 
SPIN VALVE STRUCTURES 

More physical insight can be gained by considering, an 
equivalent resistor network of the spin valve deviceJlj In 
the linear transport regime, where the measured volt- 
ages are linear functions of the applied currents, the spin 
transport for the conventional and non local geometry 
can be represented by a two terminal and four terminal 
resistor network respectively. This is shown in Fig. |^ 
for both parallel and anti-parallel configuration of the 
ferromagnetic electrodes. The resistances i?| and i?| 
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Conventional 



Non local 



Ri U 



R 



R^ 



ir 



Rt Rj, 




R 



Rj. Rt 



rt 



v- 



FIG. 2: The equivalent resistor networks of the spin valve 
device, (a) The conventional spin valve geometries in parallel 
and (b) in anti-parallel configuration, (c) The non-local spin 
valve geometry in parallel and (d) in anti-parallel configura- 
tion. 



represent the resistances of the spin up and spin down 
channels, which consists of the different spin-up and spin- 



down resistance of the ferromagnetic electrodes (R F 



T 

and the spin independent resistance R N of the nonmag- 
netic wire in between the ferromagnetic electrodes. From 
resistor model calculations we obtain: 



R\ 



= R F + R N = 



2X f 



Rf 



R N = 



w(l — ap) 



■R F H 

w 



R n 



w 



□ > 



(18) 
(19) 



where i?g = 1/aph and Rjj = arc the "square" 

resistances of the fcrromagnet and non-magnetic metal 
thin films, w and h are the width and height of the 
nonmagnetic metal strip. The resistance R — (Ajy — 
L/2)2R^/w in Fig. ||c and Fig. [2|d represents the resis- 
tance for one spin channel in the side arms of the nonmag- 
netic metal cross over a length Xn — L/2, corresponding 
to the regions IV and V of Fig. [l]b. 

Provided that Ajy 3> L the spin dependent resistance 
^ftConv k e t ween the parallel (Fig. ||a) and anti-parallel 
(Fig. |^b) resistor networks for the conventional geometry 



can be calculated using eqs. 
familiar expression:Q'Q 



L8 and 19. We obtain the 



AR 



Conv 



2(R 1+ R L )- 



(20) 



For the non local geometry and under the condition 
A at 3> L the spin dependent resistance AR between the 
parallel (Fig. ||c) and anti-parallel (Fig. |^d) resistor 
network can also be calculated. We obtain: 



AR 



m+RO' 



(21) 



Equation |T] again shows that the spin signal mea- 
sured in a non local geometry is reduced by a factor 2 
as compared to a conventional measurement. Provided 
that R F , R F <C R N we can use Eqs. |l8| and |l^ to rewrite 
Eq. |2l| into: 



AR 



2a F X F R^ 
(l-a 2 F ) 2 LwR[ 



(22) 



Using S — wh and replacing the square resistance by 
the conductivities Eq. £2] reduces to Eq. [l4|. A direct 
relation can now be obtained between the experimentally 
measured quantities AR, R^, R^ and the relevant spin 
dependent properties of the ferromagnet: 



R\ — i?t 



8ARR, 



L 4«fAf-Rq 
w (1 — a F )w 



(23) 



V. SAMPLE FABRICATION AND 
MEASUREMENT GEOMETRY 

We use permalloy Ni so Fe 2 o (Py), cobalt (Co) and 
nickel (Ni) electrodes to drive a spin polarized current 
into copper (Cu) or (Al) crossed strips. Different aspect 
ratios of the rectangular ferromagnetic injector (Fl) and 
detector strips (F2) result in different switching fields of 
the magnetization reversal process, allowing control over 
the relative magnetization configuration of Fl and F2 
(parallel/anti-parallel) by applying a magnetic field 
allel to the long axis of the ferromagnetic electrodes. t9e-3 
Two sets [F1,F2] of different sizes are used in the exper- 
iments. One set has dimensions of 2 x 0.8 /im 2 (Fl) and 
14 x 0.5 fim 2 (F2), whereas the other set has dimensions 
of 2 x 0.5 \xm 2 (Fl) and 14 x 0.15 \im 2 (F2). An example 
of a typical device is shown in Fig. 0. 

The devices are fabricated in two steps by means of 
conventional e-beam lithography with PMMA resist and 
lift-off technique. To avoid magnetic fringe fields, the 
ferromagnetic electrodes are deposited first on a ther- 
mally oxidized silicon substrate. The 40 nm thick Py 



6 




FIG. 3: a, Scanning electron microscope (SEM) picture of 
the lateral mesoscopic spin valve device with a ferromagnetic 
electrode spacing L = 500 nm. The two horizontal strips are 
the ferromagnetic electrodes Fl (Pyl) and F2 (Py2). Their 
sizes are (2x0.5 /im 2 ) and (14/imx0.15/im 2 ) respectively. An 
aluminum (Al) cross is placed in between the Py electrodes, 
which vertical arms lay on top of the Py electrodes. A total 
of 10 contacts (not all visible) are connected to the device, 
b, The conventional measurement geometry and (c) the non 
local measurement geometry. The black arrow indicates the 
direction of the applied magnetic field B in the measurements. 



electrodes are sputter deposited on a 2 nm tantalum 
(Ta) adhesion layer. The base pressure of the sputter 
system was 2 • 10~ 8 mbar vacuum, whereas the back- 
ground Ar pressure during sputtering was 1 mbar. A 
small B-ficld of 3 mT along the long axis of the Py elec- 
trodes was applied during growth in order to align the 
magneto crystalline anisotropy with the magnetic shape 
anisotropy of Py electrode. The conductivity of the Py 
film was determined to be op y — 6.6 • 10 6 fi -1 ™ -1 
and (j py = 1.2 • 10 7 ti^nr 1 at RT and 4.2 K re- 
spectively. The 40 nm thick Co (99.95 % pure) and 
30 nm thick Ni (99.98 % pure) electrodes were de- 
posited by e-gun evaporation in a 1 • 10~ 6 mbar vac- 
uum (base pressure: 2 • 10~ 7 mbar). The conductiv- 
ities of the Co and Ni films were determined to be 
a Co = 4.2 • 10 6 fl^m^ 1 and a Ni = 7.6 • 10 6 Q^m" 1 at 
RT, whereas at 4.2 K they were a Co = 7.3 • 10 6 Q- 1 m~ 1 
and aNi = 1-6 • 10 7 Q to -1 . In the second fabri- 
cation step, 50 nm thick crossed Cu (99.99 % pure) 
or Al (99.999 % pure) strips were deposited by e-gun 
evaporation in a 1 • 10~ 8 mbar vacuum (base pressure: 
2 • 10~ 9 mbar). Prior to the Cu or Al deposition, a 
few nm of Py, Co or Ni material was removed from 
the ferromagnetic electrodes by Kaufmann sputtering at 



500 Volt for 30 seconds in a 2 ■ 10~ 4 mbar Ar pressure, 
thereby removing the oxide to ensure transparent con- 
tacts. The time in between the Kaufmann sputtering 
and Cu or Al deposition was about 3 minutes. The con- 
ductivities of the Cu and Al films were determined to be 
a Cu = 3.5 • 10 7 fl^m- 1 and a Al = 3.1 ■ 10 7 ft" 1 ™" 1 at 
RT, whereas at 4.2 K they were a Cu = 7.1 • 10 7 SI -1 ™ -1 
and o M = 8.0 ■ 10 7 ft^m -1 . 

Two different measurement geometries are used to 
measure the spin valve effect in our device structure, the 
so called 'conventional' geometry (Fig j^b) and 'non-local 
geometry (Figj^c). In the conventional measurement ge- 
ometry the current is sent from contact 1 to 7 and the 
signal R = V/I is measured between contacts 4 and 9, see 
Fig. |^a. The conventional geometry suffers from a rel- 
atively large background resistance as compared to the 
spin valve resistance. Small parts of the ferromagnetic 
electrodes underneath the vertical Cu or Al wires of the 
cross are included in this background resistance, which 
can give rise to AMFo contributions and Hall effects. 
In the non-local measurement geometry the current is 
sent from contact 1 to 5 and the signal R = V/I is mea- 
sured between contacts 6 and 9. This technique is similar 
to the " potentiometric" method of Johnson used in Ref. 
|2{|[30| However, the separation of the current and voltage 
circuits, allow us to remove the AMR contribution and 
Hall cffecs of the ferromagnetic electrodes completely: 
the (magneto) resistance of the current injecting contact 
(Fl) is not relevant because any voltage drop that devel- 
ops across it will not influence the current that is sent 
through it and similarly, no current flows through the 
ferromagnetic voltage contact (F2), so its (magneto) re- 
sistance does not effect the voltage measurement. There- 
fore the only region which could possibly give rise to a 
magneto resistance signal is the nonmagnetic metal cross, 
but explicit measurements confirm that this region does 
not produce any magneto resistance. Moreover, this cen- 
tral region would be unable to produce a Hall signal in 
a non local measurement with a magnetic field perpen- 
dicular to the substrate plane as both voltage probes are 
attached longitudinal to the direction the current. 



VI. SPIN ACCUMULATION IN PY/CU/PY 
SPIN VALVES 

The measurements were performed by standard ac- 
lock-in-techniques, using current magnitudes of 100 \iA 
to 1 mA. Typical spin valve signals of two samples MSV1 
and MSV2 (of the same batch) with a Py electrode spac- 
ing of L — 250 nm are shown in the Figs. UJI and 0. 
They are both measured in a non local measurement ge- 
ometry and conventional measurement geometry. Sam- 
ple MSV1, data shown in Fig. |] and ||, had a current 
injector Pyl electrode of size 2 x 0.5 fj,m 2 , whereas the 
detector electrode Py2 had a size of 14 x 0.15 nm 2 . Sam- 
ple MSV2, data in shown Fig. ^], had wider Py electrodes 
of 2 x 0.8 [im 2 and 14 x 0.5^m 2 . The first set of (nar- 
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rower) Py electrodes [Pyl,Py2] had a more ideal switch- 
ing behavior and had three times larger switching fields 
as compared to the second Set [Pyl,Py2]. We note that a 
discussion of the magnetic behavior of the Py electrodes 
and contacts has been given in Ref. [l9| To control the 
parallel/anti-parallel magnetization configuration of the 
Py electrodes a magnetic field is applied parallel to the 
long (easy) axis of the Py electrodes. 



A. Non local spin valve geometry 

Figure ^a and [|b show typical data in the non local 
measurement geometry taken at 4.2 K and RT for sample 
MSV1 with a 250 nm Py electrode spacing. Sweeping the 
magnetic field from negative to positive field, an increase 
in the resistance is observed, when the magnetization of 
Pyl flips at 9 mT, resulting in an anti-parallel magneti- 
zation configuration. The rise in resistance is due to the 
spin accumulation or equivalently an excess spin density 
present in the Cu metal. When the magnetization of Py2 
flips at 47 mT (T = 4.2 K) and 38 mT (RT), the magneti- 
zations are parallel again, but now point in the opposite 
direction. The magnitude of the measured background 
resistance, around 30 mVt at T = 4.2 K and 120 mfi at 
RT, depends on the geometrical shape of the Cu cross 
and is typically a fraction of the Cu square resistance. 

Figure [|c and |]d show the " memory effect" . Coming 
from high positive B field, the sweep direction of the B 
field is reversed after Pyl has switched, but Py2 has not. 
At the moment of reversing the sweep direction, the mag- 
netic configuration of Pyl and Py2 is anti-parallel, and 
accordingly a higher resistance is measured. When the 
B-ficld is swept back to its original high positive value, 
the resistance remains at its increased level until Pyl 
switches back at a positive field of 9 mT. At zero B field 
the resistance can therefore have two distinct values, de- 
pending on the history of the Py electrodes. 



B. Conventional spin valve geometry 

The top curve in Fig. |5| shows the magneto resistance 
behavior of sample MSV1 in the conventional measure- 
ment geometry. A small AMR contribution (dip in curve) 
of the Pyl electrode around —9 mT and a small Hall 
signal caused by the Py2 electrode can be observed in 
the negative sweep direction. Because a small part of 
the Py electrodes underneath the Cu wire is measured 
in this geometry, (local) changes in the magnetization 
at the-Ey/Cu contact area can produce an AMR or Hall 
signal.Ej In the positive sweep direction a dip is no longer 
observed indicating that the magnetization reversal of 
the Pyl electrode is not the same for the two sweep di- 
rections. However, in the magnetic field range in be- 
tween the two switching fields, we do observe a resistance 
'plateau' from 10 mT up to a field of 45 mT. 
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FIG. 4: The spin valve effect at T = 4.2 K (a) and RT (b) 
in the non-local geometry for a Py/Cu/Py spin valve device 
(sample MSV1) with 250 nm Py electrode spacing. An in- 
crease in resistance is observed, when the magnetization con- 
figuration is changed from parallel to anti-parallel. The solid 
(dashed) lines correspond to the negative (positive) sweep di- 
rection. (c),(d) illustrate the "memory effect". For clarity the 
(c) and (d) are off set downwards. Note that the vertical scale 
of (a) is different from (b),(c) and (d). The sizes of the Pyl 
and Py2 electrodes are 2 x 0.5 /xm 2 and 14/im x 0.15/xm 2 



The magnitude of the spin valve effect measured in 
the conventional geometry is about 4.1 mVL at T = 4.2. 
This is about 2.5 times bigger than the magnitude of the 
spin signal measured in a 'non-local geometry (1.6 mfi). 
Note that the factor 2.5 is deviating from the factor 2 as 
predicted by the eq. 17. This is attributed to deviations 
from our 1-dimensional model, which can be expected 
for the samples with the shortest Py electrode spacing 
L = 250 nm, as the presence of the Cu side arms for 



these samples, see Fig. t 
The top curve in Fig. 



are most felt, 
j shows the magneto resistance 
behavior in the conventional measurement geometry for 
the sample MSV2. Here a change of the resistance is 
already observed before the field has reached zero in a 
positive field sweep, whereas the negative field sweep is 
very asymmetrical compared to the positive field sweep. 
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FIG. 5: The spin valve effect of sample MSV1 in a con- 
ventional measurement geometry (top curve) at T — 4.2 K 
and non-local measurement geometry (bottom curve), with 
a Py electrode spacing L = 250 nm. The sizes of the Py 
electrodes are 2 x 0.5 /im (Pyl) and 14 x 0.15 /im 2 (Py2). 
The solid(dotted) curve corresponds with a negative (posi- 
tive) sweep direction of the B-field. 



FIG. 6: The spin valve effect of sample MSV2 in a con- 
ventional measurement geometry (top curve) at T = 4.2 K 
and non-local measurement geometry (bottom curve), with 
a Py electrode spacing L = 250 nm. The sizes of the Py 
electrodes are 2 x 0.8 fxm (Pyl) and 14 x 0.5 /im 2 (Py2). 
The solid(dotted) curve corresponds with a negative (posi- 
tive) sweep direction of the B-field. 



This is attributed to the formation of a multi-domain 
structure in the 2 x 0.8 /im 2 (Pyl) electrode, causing a 
large AMR (w 10 mO) signal at the Py/Cu contact area 
of the Pyl electrode. 

However, in a non local measurement geometry, the 
" contact" magneto resistance contribution of the Py elec- 
trodes can be removed and a clear spin valve signal is ob- 
served with a similar magnitude as sample MS VI. This 
is shown in the bottom curve of Fig. ||. Note also that 
the larger widths and aspect ratio of the Py electrodes 
in sample MSV2 result in three times smaller switching 
fields as compared to sample MS VI. 



C. Dependence on Py electrode spacing 

A reduction of the magnitude of spin signal AR is ob- 
served with increased electrode spacing L , as shown in 
Fig. H By fitting the data to Eq. |l| we have obtained 
the spin relaxation length Xn in the Cu wire. From the 
best fits we find a value of 1 fj,m at T = 4.2 K, and 
350 nm at RT. These values are compatible with those 
reported in literature, where 450 tun is obtained for Cu 
in GMR measurements at 4.2 K.o However a detailed 
discussion on the obtained spin relaxation lengths and 
corresponding spin relaxation times will be given in Sec. 

ixl 



In principle the fits of Fig. |7] also yield the spin po- 
larization ap and the spin flip length Xf of the Py 
electrodes. However, the values of and Xf can- 
not be determined separately, as in the relevant limit 
(M >> 1) which applies to the Py/Cu/Py experiments 
(12 < M < 26), the spin signal AR is proportional to the 
product o>fXf as is shown by Eq. [l4|. From the fits we 
find that ofXf = 1-2 nm at 4 , 2 K a nd o.fXf = 0.5 nm 
at RT. Taking, from liter atureE3E30, a spin flip length in 
the Py electrode of Xp — 5 nm (at 4.2 K), a bulk current 
polarization of ~ 20 % in the Py electrodes at T = 4.2 
K is obtained: ap = 0.2. We note however that the 
injected spin polarized current from the Py electrode is 
partially shunted by the Cu wire lying on top of the Py 
electrode. When taken into account we estimate that it 
could increase the value upXp by a factor 2 to 3. 

It is also possible to calculate the polarization of the 
current at the Py/Cu interface. For a sample with a 
Py electrode spacing of L = 250 nm at T = 4.2 K and 
using Eq. [l5| we find: P ~ 0.02, a factor 10 lower than 
the bulk polarization ap of the Py electrodes. From the 
resistor model we can see why the current polarization 
at the Py/Cu interface is reduced. For this we need to 
calculate the magnitude of the spin dependent resistance 
difference. Using Eq. |2^ and L — 250 nm, AR = 1.6 mfi, 
= 0.3 fi and w = 100 nm ( at T = 4.2 K) we find: 
Rl — i2| w 100 mfi. From the right hand side term of Eq. 
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FIG. 7: Dependence of the magnitude of the spin signal AR 
on the Py electrode distance L, measured on Py/Cu/Py sam- 
ples in the non local geometry. The solid squares represent 
data taken at T = 4.2 K, the solid circles represent data 
taken at RT. The solid lines represent the best fits based on 
equation ll2L 



and using Rf, 



2 SI we can check that this indeed 



corresponds with the value of apXp rj 1.2 nm, as was 
also obtained from the fit in Fig. ^. From Eqs. 
and using Xp = 5 nm, ctF — 0.2 (at 4.2 K) we obtain the 
spin up and spin down resistance of the Py ferromagnet: 



R?y = 



2Xi 



R 



Py _ 



w(l + ap ) 

2X F 
w(l — ap) 



■Rf 



Rf 



160 mf] 



260 mfi 



(24) 
(25) 



This shows that the total resistance experienced over 
a length Xf + Xn by the spin up and spin down currents 
is indeed dominated by the spin independent resistance 



R + R = Xn2R^/w ~ 6 CI. Here we have used that 
Aw = 1 fim at T = 4.2 K and w — 100 nm. This leads to 
a interface polarization of P w (i? J .-i? T )/(i? Ar + R) m 2% 
at the Py/Cu interface. 

All though the role of interface resistance for spin in- 
jection will be described in the next section, we note here 
that the small difference R± — i2j « 100 mCl responsible 
for a spin valve signal of AR = 1.6 mCl could possibly 
also result from an interface resistance at the Py/Cu in- 
terface. Commonly rep»eted-*csistjvitics of 5 • 10~ 16 Clm 2 
for the Py/Cu interfaceOt3E3c3e3 and a contact area of 
S = 1 ■ 10~ 14 m 2 (i.e. Rint — 50 mCl) would yield a 
realistic interface polarization of 77 = 0.5 for the Py/Cu 
interface, using eq. ||(]. However, the specific details of 
the spin injection mechanism (interface, bulk or a com- 
bination) do not alter the conclusion that the total spin 
dependent resistance R± — Rf « 100 mCl is dominated 
by the spin independent resistance of the Cu strip over 
a spin relaxation length and hence leads to a consider- 
able reduction of the spin valve signal, as was pointed 
out above. 



D. Comparison with Johnson Spin Transistor 

The magnitude of the spin signals in the Py/Cu/Py 
samples, when scaled to the cross sections utilized in 
the Au thin film devices of Ref. p9||3C| (the "Johnson 
spin transistor"), are more than lO^imes smaller than 
obtained in that previous work. However, in that ear- 
lier work it was necessary to invoke a spin polarization 
exceeding 100i%|-to explain the results in terms of spin 
accumulation.E3E!l This contrasts our results, which yield 
a spin polarization P of the current injected in the Cu 
wire at the Py/Cu interface of about 1 — 2%. 

In Refs. ||^,^,^8| and ^ Johnson postulates that spin 
injection is mediated by interfacial transport, because the 
interface resistances i?™*, RY 1 * would dominate the to- 
tal resistance in both spin- up and spin-down channels: 



lit 



> 



Rf 



R 



N 



R and R 



> 



Rf 



R + R respec- 



tively. In this limit spin injection would be characterized 
by the interfacial spin injection parameter defined as: 



V = 



TDint 



r>int 



(26) 



and Johnson derives 
spin accumulation signal 



AR = 



following expression for the 



2 V 2 X 



N 



un SL 



(27) 



Applying Eq. |27] Johnson calculates an expected spin 
signal of AR = 1.9 £1 for our device with the shortest Py 
electrode spacing L = 250 nm, using S = 5 • 10~ 15 m 2 , 
a Cu = 7.1 • 10 7 f^ 1 ™ -1 , r] = 0.4 and X N = 1.0 fimH 

However a polarization of the current at the Py/Cu 
interface of 7/ = 40 % would require spin dependent in- 



terface resistances of R 



16 tt and RV™ = 37 ft, using 



Eqs. [2l] and replacing Eqs. [l8] and [l 7 ^ by 



Rf 



R 



N 



R 1 = R\ 



R\ +R 



N 



(28) 
(29) 



where the spin dependent interface resistances i?™* and 
have simply been added up to bulk spin dependent 

resistances Rf and Rf because the spin polarization 77 
as well as the bulk spin polarization ap are foij 
positive (up > and 77 > 0) for Py and Cui 
The values Rf 1 = 16 fi and Rf 1 = 37 CI yield a total 
single interface resistance Ri n t = 11 CI or equivalently, 
a interface resistivity of 1 • 10~ 13 Ctm 2 . This is more 
than a 100 times larger then the upper limit 0.1 SI or 
equivalently a contact resistivity of 1 • 10~ 15 dm 2 that 
we are able to determine from our experiment, see Figs. 
I and |. 

The above arguments also apply for the experiment 
of Refs. ^|[50| where a gold layer is sandwiched in be- 
tween two Py layers. There is no physical reason why 
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there should exist an interface resistivity larger than 
1 • 10~ 13 Vim 2 between the Au and Py or Co layers in 
the experiment of Ref. ^9| which can explain an inter- 
face current polarization of r\ — 0.4 or more. Equation 
27 can therefore not be applied to the experiment of Ref. 
29, because it does not include the (fast) spin relaxation 
reservoirs of the ferromagnetic injector and detector con- 
tacts, which dominate the total spin relaxation in the 
case of transparent contacts, as was already pointed out 
in Refs. 

In view of this, given the unexplained discrepancies 
(77 > 3) of the earlier work in Ref. ^9|,^0[ and the more 
consistent values obtained in the recent work, it is our 
opinion that the results of Refs. 29 3^ cannot be recon- 
ciled with spin injection and spin accumulation. 



VII. SPIN ACCUMULATION IN PY/AL/PY 
SPIN VALVES 



Here we will describe spin injection experiments using 
permalloy Ni$oFe2o (Py) strips as ferromagnetic elec- 
trodes to drive a spin polarized current via transparent 
contacts into aluminum (Al) crossed strips, see Fig. |^. 
Similar current polarizations and spin relaxation lengths 
for P y a nd Al are obtained as in the previous section 
(Sec. 0). 



A. Spin valve measurements 

Figure || shows a typical spin valve signal of a 
Py/Al/Py sample with a Py separation spacing of L = 
250 nm and Py electrodes of sizes 2 x 0.8 fim 2 and 
14 x 0.5 /iTO 2 . 

The top curve in Fig. || shows the magneto resis- 
tance behavior in the conventional measurement geome- 
try. Again the magneto resistance signals of the Py con- 
tacts are dominating in this geometry, reaching a maxi- 
mal amplitude of about 6 mfL Note that the two resis- 
tance values at high positive and negative field differ by a 
value of about 0.3 mfi, which is attributed to a local hall 
effect caused by the 14 x 0.5 /dm 2 Py electrode. The bot- 
tom curves in Fig. ^ show magnetic field sweeps in the 
non local measurement geometry, which clearly shows a 
spin valve signal having removed all the spurious contact 
magneto resistance effects. The magnitude of the spin 
valve signal measured is 0.18 mfl at 4.2 K. 



B. Dependence on Py electrode spacing 

A reduction of the magnitude of spin signal AR of the 
Py/Al/Py samples is observed with increased electrode 
spacing L, as shown in Fig. ||. However, for the T = 4.2 
K data this dependence is not monotonic. The spin valve 
devices with small L — 250 nm and L — 500 nm show a 
smaller spin valve signal than the device with L = 1 fim. 
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FIG. 8: The spin valve effect of a Py/Al/Py sample using 
a conventional measurement geometry (CONV, top curve) at 
T = 4.2 K and non-local measurement geometry (NL, bottom 
curve), with a Py electrode spacing L = 250 nm. The sizes 
of the Py electrodes are 2 x 0.8 /im 2 (Pyl) and 14 x 0.5 /Ltm 2 
(Py2). The solid(dotted) curve corresponds with a negative 
(positive) sweep direction of the B-field. 



We note that all the devices shown in Fig. |9| are from the 
same (processing) batch. However, the granular struc- 
ture of the Al film with a grain size in the order of the 
width of the Al strip causes fluctuations in the resistance 
of the Al strip in between the Py electrodes. The sample 
with L = 250 nm and L = 500 nm indeed show a higher 
resistance than expected when measured in the conven- 
tional geometry at T = 4.2 K. This irregular behavior of 
the resistance due to grains is not observed at Rt due to 
the additional presence of phonon scattering. From the 
best fits to Eq. [l2| we find a spin relaxation length 
in Al of 1.2 nm at T = 4.2 K and 600 nm at RT. Note 
that the spin flip lengths are about 2 times larger than 
reported in Ref. The reason for this increase is the 
higher conductivity of the Al in these samples, caused by 
a lower background pressure of 1 • 10 -8 mbar during evap- 
oration as compared to a background pressure of 1 ■ 10~ 6 
used in Ref. 
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The fits of Fig. || also yield the spin polarization cv^ 
and the spin flip length Xp of the Py electrodes. We 
find apXp — 1.2 nm at 4.2 K and apXp = 0.5 nm 
at RT, in agreement with the Py/Cu/Py spin valve data 
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FIG. 9: Dependence of the magnitude of the spin signal AR 
on the Py electrode distance L, measured in the non local 
geometry for Py/Al/Py spin valves. The solid squares repre- 
sent data taken at T = 4.2 K, the solid circles represent data 
taken at RT. The solid lines represent the best fits based on 
equation Fl2l 



of Sec. [v| Note that for the Py/Al/Py spin valve also 
applies that M >> 1 and thus the spin signal AR 
is proportional to the product cxfXf (25 < M < 32). 
Using Eq. |l5|, a polarization P for the Py/Al/Py sample 
with the smallest Py electrode spacing of L = 250 nm 
at T = 4.2 K is found to be only 3%: P = 0.03. 



The magneto resistance traces of Fig. |10| a indicate a 
clear switching of the magnetization at 20 mT of the 
14 x 0.5 tim 2 Co2 electrode and is attributed to a local 
Hall effect produced at the Co/Cu contact area of the 
Co2 electrode. 
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VIII. SPIN INJECTION USING CO AND NI 
FERROMAGNETIC ELECTRODES 



From Eq. 



it can be seen that the magnitude of 



the spin dependent resistance Ri — i?| is sensitive to the 
properties ap, Xf and op of the ferromagnet. As R\ —R\ 
enters squared in the spin valve signal AR, see Eq. |2l|, an 
increase of Xf with a factor 10 would increase AR with 
a factor 100. We have therefore tried cobalt (Co) and 
nickel (Ni) as ferromagnetic spin injectors and detectors 
to increase the magnitude of the spin valve signal, as 
larger spiB,,£elaxation lengths can be expected for these 
materials, tm 



A. Spin accumulation in Co/Cu/Co spin valves 

Figure |Io| a shows a " contact" magneto resistance trace 
and magnetic switching behavior at RT of a 14 x 0.5 /im 2 
(Co2) electrode of a Co/Cu/Co spin valve device with a 
Co electrode spacing of 250 nm and Co electrodes of sizes 
2 x 0.8 /im 2 and 14 x 0.5 fim 2 . The "contact" magneto 
resistance is measured by sending current from contact 5 
to 7 and measuring the voltage between contacts 6 and 
9 (see Fig. |3|). Note that in this geometry the measured 
voltage is not sensitive to a spin valve signal as only one 
Co electrode is used in the measurement configuration. 



FIG. 10: a, "Contact" magneto resistance trace of the Co2 
electrode with size 14 x 0.5 fj,m 2 . The Hall signal indicates 
an abrupt magnetization switching of the Co2 electrode, b, 
The spin valve effect at RT in a Co/Cu/Co device with a Co 
electrode spacing L — 250 nm, using the non local measure- 
ment geometry. The solid (dotted) curve corresponds with a 
negative (positive) sweep direction of the B-field. 



Figure |10jb shows the spin valve effect at RT for a 
Co/Cu/Co spin valve device. The magnitude of the spin 
dependent resistance AR — 0.25 mfl is slightly smaller 
than in the Py/Cu/Py spin valve device. At T = 4.2 
K the signal increases to AR = 0.8 mVt. Using Eq. |l2| 
and the values of cat, Atv for Cu and <jp for Co (see Sec. 
0), we obtain olf\f = 0.3 at RT and cxf^f — 0.7 at 
T = 4.2 K. These obtained values are much smaller than 
reported for Co in GMR, r^qyrjrrjfvnts where a? w 0.5 



and Xf — 10 — 60 nm 
be discussed in Sec 



iff 



IXC 



This discrepancy will 



B. Spin accumulation in Ni/Cu/Ni spin valves 

In Fig. [Tl] a and b two "contact" magneto resistance 
traces of a Ni electrode (Nil) with size 2 x 0.5 /xm 2 (top 
curve) and a Ni electrode (Ni2) with size 14 x 0.15 /iro 2 
(middle curve) are shown of a Ni/Cu/Ni spin valve de- 
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FIG. 11: (a) "Contact" magneto resistance trace (see text) 
of the Nil electrode with size 2 x 0.5 /im 2 . (b) "Contact" 
magneto resistance trace of the Ni2 electrode with size 14 x 
0.15 fj,m 2 . (c) The spin valve effect of a Ni/Cu/Ni device 
at RT with a Ni electrode spacing of L = 500 nm, using a 
non local measurement geometry. The solid(dotted) curve 
corresponds with a negative (positive) sweep direction of the 
B-field. 



IX. SPIN RELAXATION TIMES OF 
CONDUCTION ELECTRONS IN METALS 

In this section we will analyze our obtained spin re- 
laxation times r s f in CuandAl from the spin injection 
experiments in Sees. VI, VII VIII and compare them with 
theorjt-and previously reported values from CPP-GMRE3, 
CESRcH, weak localization!^ and superconducting tun- 
neling experimental The obtained spin polarization and 
spin relaxation lengths in Py, Co and Ni will be compared 
with reported values from CPP-GMR experiments. 

In CESR experiments the measured electron spin 
transverse relaxation time Ti is proportional to the 
width_of the absorption peak at the resonance frequency. 
YafetO showed that in metals T2 is equal to the longitudi- 
nal spin relaxation time or spin flip time T± (T\ = r s f). In 
weak localization and superconducting tunneling experi- 
ments the spin orbit scattering time r s o . is determined, 
with Tg. . being defined similarly in both experiments.t2l 
Spin orbit interaction in weak localization experiments 
is responsible for destructive interference when electrons 
are scattered at (nonmagnetic) impuritiestS, whereas 
in the superconducting tunneling experiments it mixes 
up the spin-up and spin-down quasi-particle density of 
states, wlign_they are Zeeman-split by an applied mag- 
netic fieldE'Eil We make the identification r s o . = r s /. 



Discussion of electron spin relaxation in 
nonmagnetic metals 



vice with a Ni electrode spacing of 500 nm. For the Nil 
contact current is send from contact 1 to 5 and the volt- 
age is measured from contact 4 to 6 (see Fig. ||). For the 
Ni2 contact current is send from contact 5 to 7 and the 
voltage is measured from contact 6 to 9. In the magnetic 
field sweeps of Figs. [il] a and [ll| b a large range can be 
observed where the magnetization configuration of the 
Ni electrodes are anti-parallel. However no spin valve 
signal could be detected within experimental accuracy in 
the non local measurement geometry at RT as well as at 
T = 4.2 K, as is shown in Fig. [TT| c (RT). An upper bound 
on the spin valve signal is found to be AR < 20 fi£l at 
RT as well as at T = A.2K. Using Eq. 12 and the values 
of cat, Aat for Cu and up for Ni (see Sec. |v|), we ob- 
tain apXp < 0.3 at RT as well as at T = 4.2 K. These 
obtained values are also much smaller than reported for 
Ni in GMR experiments, whstfr-etp w 0.2 and using the 
expected Xp(calc) = 15 nm.EJ'E3 This discrepancy will 



be discussed in Sec. IX C 



We note that the magnetic field in the measurements 
of Fig. [n] is applied perpendicular to the long axis of 
the Ni electrodes, showing a more pronounced magnetic 
switching behavior than an applied magnetic field along 
the long axis of the Ni electrodes. 



The fact that a spin can be flipped implies that there 
is some mechanism which allows the electron spin to in- 
teract with its environment. In the absence of magnetic 
impurities in the nonmagnetic metal, the dominant mech- 
anism that provides for this interaction is the spin-orbit 
interaction, as was argued by Elliot and Yafet.HS When 
included in the band structure calculation the result of 
the spin-orbit interaction is that the Bloch eigenfunc- 
tions become linear combinations of spin-up and spin- 
down states, mixing some spin-down character into the 
predominantly spin-up states and vice versaS Using a 
perturbative approach Elliot showed that a relation can 
be obtained between the elastic scattering time (r e ), the 
spin relaxation time (t s /) and the spin orbit interaction 
strength defined as {X/ AE) 2 : 



T e , A 2 

T s f AE 



(30) 



where A is the atomic spin-orbit coupling constant for 
a specific energy band and AE is the energy separation 
from the considered (conduction) band to the nearest 
band which is coupled via the atomic spin orbit inter- 
action constant. Yafet has shown that Eq. is tem- 
perature independents. Therefore the temperature de- 
pendence of (t s /) _1 scales with the temperature behav- 
ior of the resistivity being proportional to (t^ 1 ). For 
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many clean metals the temperature behavior of the re- 
sistivity is dominated by the electron-phonon scattering 
and can to a good approximation be described by the 
Bloch-Griineisen relationEfl: (t s /) _1 ~ T 5 at tempera- 
tures below the Debye temperature To and ~ T 
above To. Usin g data from CESR experiments, Monod 
and BeuneuE3C3 showed that (r,/)" 1 follows the Bloch- 
Griineisen relation for monovalent alkali and noble met- 
als. In Fig. [l^ their results are replotted for Cu and 
Al, using the. revised scaling as applied by Fabian and 
Das SarmaS In addition we have plotted the obtained 
data points for Cu and Al at T /Tp ~ 1 from our spin in- 
jection experiment by calculating r s / from A s / (see Sec. 
KB below) and using the calculated spin orbit strength 
parameters from Ref. || (X/AE) 2 = 2.16 • 1(T 2 for Cu 
and {X/AE) 2 = 3 • lO^for Al. 




FIG. 12: The (revised) Bloch-Griineisen plotEI. The quantity 
C ■ {rl h s Y x is plotted versus the reduced temperature T/Td 



on logarithmic scales 
' -l 



C represents a constant which links 



( r f/ ) 1 to the (original) plotted width of a CESR resonance 
peak, normalized by the spin orbit strength (X/AE) 2 and 
the resistivity p D at T = T D : C = (-/(X/AEfpo)' 1 . Here 
(r^f ) _1 is the phonon induced spin relaxation rate, 7 is the 
Larmor frequency and Tb is the Debye temperature. We used 
p D = 1.5-10" 8 andr3jpri= 315 K for Cu and p D = 3.3-10" 8 and 
Td = 390 K for A1.E3&I The dashed line represents the general 
Bloch-Griineisen curve. The open squares represent Al data 
taken from CESR and the JS spin injection experiment (Refs. 
|lf],[3^) . The open circles represent Cu data taken from CESR 
experiments (Refs. ^,^,^). The solid square (Al) and circle 
(Cu) are values from the spin injection experiments described 
here and Refs. lilE! 



12 it can be seen that for Cu the Bloch- 
including the newly 



From Fig. 
Griineisen relation is well obeyed 
added point deduced from our spin injection experiments 
at RT (T/T D = 0.9). For Al however the previously ob- 
tained data points as well as the newly added point from 
the injection experiments at RT (T/Td = 0.75) are de- 
viating from the general curve, being about two orders 
of magnitude larger than the calculated values based on 
Eq. |3(]. We note that we cannot extract data for the 
Bloch-Griineisen plot shown in Fig. [l2] from our spin 
injection experiments at T = 4.2 K, because the impu- 



rity (surface) scattering rate is completely dominating 
the phonon contribution at T — 4.2 K. 

Fabian and Das Sarma have resolved this discrepancy 
for Al by pointing out that there can exist so called 'spin- 
hot-spots' at the Fermi surface of poly valent metals(likc 
Al) . Performing an ab initio pseudopotential band struc- 
ture calculation of Al they showed that the spin flip 
contribution of these (small) spin-hot-spot areas on the 
(large) Fermi surface dominate the total spin flip scatter- 
ing rate r s / , making it factaLlQO—faster than expected 
from the Elliot- Yafet relationOO'B A simplified reason- 
ing for the occurrence of these spin-hot-spots is that in 
poly valent metals the Fermi surface can cross the first 
Brillioun zone making the energy separation AE in Eq. 
|30| between the (conduction) band to the spin orbit cou- 
pled band much smaller at these (local) crossings and 
hence result in a larger spin orbit strength (X/AE) 2 . 
Our newly added data point shows that the under es- 
timation of the spin orbit strength also holds for Al at 
RT (T/T D = 0.75), as can be seen in Fig.|f. However it 
is in excellent agreement with the theoretical predicted 
value by Fabiaro as will show below. 



B. Qualitative analysis of the spin relaxation time 
r s f in Cu and Al 



Comparing the conductivities and spin relaxation 
lengths at RT and T = 4.2 K we can obtain the impu- 
rity and phonon scattering rate and their associated spin 
relaxation rates. Therefore we can define an impurity 
spin relaxation ratio a lmp — T %mp /t 1 ™ p and an inelastic 



•sf 

(phonon) scattering ratio a ph — r ph /t v ^ . Here (t^p)- 1 
and (t p ' 1 )~ 1 are the impurity and phonon scattering rate 
and {r i J? p )~ l and (t^ 1 ) -1 are the impurity and phonon 
induced spin relaxation rate. From the measured conduc- 
tivity at T — 4.2 K and Eq. || we can determine r lmp . 
Using the Mathiessen rule (t^- 1 = (t™p)- 1 + (t^)' 1 
and the RT conductivity we can determine t v . We note 
that the surface scattering in our samples is dominat- 
ing the impurity scattering, as the mean free paths of 
l e » 60 nm for both Al and Cu at T = 4.2 K are larger 
than their film thicknesses (50 nm). In the calculation 
we use the free electron values N(Ep)c u — 1-8 • 10HE 
states/eV/m 3 and v F (Cu) = 1.57 • 10 6 m/s for CuO 
and we use N(E F ) A i = 2.4 v40 28 states/eV/m 3 and 
v F (Al) = 1.55 • 10 6 m/s for AI.O 

The obtained parameters for Cu and Al ( T ]/ lP j' s / 
a imp , a ph ) are tabulated in tabel |. From tabel | we see 
that rfj and a p h for Al at RT are in good agreement with 
the theoretical values as predicted inJdie bandstructure 
calculation by Fabian and Das Sarma.O They are also in 
agreement with the results obtained from CESR exper- 
iments and the earlier JS spin injection experiments at 
temperatures below 90 K. Note that the spin relaxation 
times are 2 orders of magnitude larger in those earlier 
experiments due the use of extremely clean samples with 
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TABLE I: Comparison of spin relaxation times between different experiments. t % J^ v [ps] is the impurity induced spin relaxation 

time at low temperatures T < 4.2 K due to surface scattering, dislocations or grain boundaries, rf* [ps] is the phonon induced 
spin relaxation time at elevated temperatures due to inelastic phonon scattering. For definition of a,i mp and a p h see text. 

"For T=293 K 
'For T=45 K 

Tor a temperature range T=[1..90] K 
d For a temperature range T=[1..60] K 



electron mean free paths of a few tens of micrometers. 
Also for Cu we see that rfj and a p h at RT are in good 
agreement with the results obtained from CESR experi- 
ments at temperatures below 60 K. 

The impurity scattering ratio a mp shows a much big- 
ger spread in values for both Al and Cu. We specu- 
late that this due to the different origin of the impurities 
in the samples used for the various measurement tech- 
niques. For the CESR experiments the impurity scat- 
tering is caused by dislocations, whereas for our exper- 
iment it is mainly due to surface scattering. In weak- 
localization and superconducting tunnneling experiments 
the used films are even thinner then our films (few nm), 
causing additional scattering from grain boundaries in 
addition to the surface scattering. 

We note that for thin films we use it is not possible 
to realize mean free paths of the order of micrometers as 
they will always be limited by surface scattering. How- 
ever the sensitivity of the CESR technique does not allow 
measurements of below typically 1 ns, whereas the 
SQUID detection technique used in Ref. 16 does not op- 
erate at RT. Therefore spin injection into thin films is 
rather complementary to the CESR techniques and the 
JS spin injection experiments in the determining r^ 1 in 
the temperature range from liquid Helium to RT. Also, 
the fact that about half of the momentum scattering pro- 
cesses at RT is due to phonon scattering implies that the 
present obtained results on the spin relaxation lengths in 



Al and Cu can be maximally improved by a factor of 2 
at RT. 



C. Spin injection efficiency of Py, Co and Ni 
ferromagnets 

In addition to the spin-orbit spin scattering in metallic 
ferromagnets, as described above for nonmagnetic met- 
als, there is spin flip scattering by magnons £3 In cobalt 
magnons are nearly absent at low temperatures and only 
start to compete with the spin orbit spin-flip scattering 
at temperatures higher than T = 100 K£j The spin flip 
scattering by magnons has two effects. It will simply 
add to the spin orbit spin flip scattering rate which re- 
duces the spin relaxation length Xp of the ferromagnet 
at higher temperatures. Secondly, it will lower the bulk 
current polarization of the ferromagnet «f by changing 
<7f and <7| and in addition by giving rise to a " spin mixing 
rate" which equalizpsiie spin-up and spin-down currents 
in the ferromagnet J 38 !! 39 ! The presence of spin flip scatter- 
ing by magnons can therefore lower of as well as Xf at 
RT. 

At low temperatures (T < lOOi-T) and in absence 
of magnetic impurities an upper estimate can be given 
for the expected spin relaxation length in Co and Ni 
due to the spin orbit spin flip scattering only: = 
\J Dr e /a, where a is taken from spin flip scattering cross- 
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sections determined by CESR experimeiitsEiKi and re- 
cently from magneto-optic experimentsB ap e = 1.1 • 
1CT 2 , a Nl = 1.5 • 1CT 2 and a Co = 4.2 • 10~ 2 . Using a 
free electron model, the spin relaxation length Xp y for 



Py with ap y 
aco = 1-7 
way in Ref. 



10 6 fi" 



= 8.1 

io 7 n^m- 1 

]5t X# y (calc) 



and Xco for Co with 



have been estimated in this 
w 9 em and Ac (ca^c) w 36 
nm at T = 4.2 K.E3 Note that Ap scales linearly with r e 
and thus the conductivity of the ferromagnetic metal. 
In this respect the reported value of \ Co — 59 nm 
in Ref.|5^ is quite remarkable, because the conductiv- 
ity of the Co metal (<r co = 6.4 • 10 6 ) used in Ref.|52| is 
about 3 times smaller then used to calculate \c {calc) 
in Ref.|35| (a co = 1.74 • 10 7 ), which makes the expected 
Xco(calc) = 13 nm in the experiment of RefJE^. For Ni 
we derive an estimate of Xp using a free electron den- 
sity of 5.4 • 10 28 m" 3 . With a Nl = 1.6 ■ 10 7 rr 1 ™" 1 
and QNi = 1.5 • 10~ 2 we calculate: \Ni(calc) = 15 nm at 
T = 4.2 K. 

Because M > 10 for all our spin valve samples, we 
cannot separately determine ap and A p from the magni- 
tude of the spin valve signal AR. In table [l| we therefore 
give the "spin injection efficiency" olfXf together with 
reported values from GMR experiments. We note that 
our thin film conductivities for Py, Co and Ni are within a 
factor 2 of the reported values in the GMR experiments. 
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TABLE II: Spin injection efficiencies cvfAf in nm for three 
different ferromagnetic metals. The data is deduced from the 
mesoscopic spin valve (MSV) experiments with transparent 
contacts in a non local geometry using Cu as nonmagnetic 
metal and compared with results from GMR experiments. 




0.2) and using Ajv;(caZe) = 15 nm 



Table || shows that our obtained spin injection effi- 
ciency of the Py ferromagnet ap y Xp y is in quantative 
agreement with the values reported in GMR experiments 
(ap y = 0.7, Xp y = 5 nm), taking into account that our 
obtained ap y Xp y represents a minimal value due to a 
partially shunting of the injected current by the Cu wire 
on top of the Py electrodes. The reduction of ap y Xp y at 
RT beyond the ratio 1.8 of the Py conductivies at T = 4.2 
K and RT could be attributed to magnons lowering ap 
at RT. 

For the Co and the Ni ferromagnets we observe much 
smaller spin injection efficiencies cxfXf, being more than 
1 order of magnitude smaller than values of cxfXf ob- 
tained in GMR experiments. So the question is, what is 
causing this rather large reduction of the the spin valve 



signal? 

First we discuss the possible influences of an exist- 
ing interface resistance at the Co/Cu and Ni/Cu inter- 
faces. From the resistance measured in a conventional 
geometry we are able to determine an upper estimate of 
the (diffusive) interface resistances. For the Co/Cu/Co 
spin valve of Fig. [io] we find an upper limit for a sin- 
gle Co/Cu interface of 0.4 fi, whereas for the Ni/Cu/Ni 
spin valve of Fig. [ll] we find for a single Ni/Cu inter- 
face 0.6 f2. We note that the associated interface re- 
sistivity (« 5 • 10~ 15 ilm 2 ) values are about 5 times 
larger thai^-cakjttlated for Co/Cu (specular or diffusive) 
interfacesuO'E3'E3 and also 5 tirocs_Larger than values ob- 
tained from GMR experimentsrllEil In case these Co/Cu 
and Ni/Cu interface resistances are spin dependent, the 
spin signal would be (largely) increased as the sign of 
the bulk and interface spin asy mmetric] s-pf Co, Ni and 
Cu are found both to be positiveE3E30E3 (a^ > and 
rj > 0). However this is clearly not observed. In the op- 
posite case of spin independent interface resistances, the 
interface resistance for each spin channel (« 1 Q) will not 
reduce the measured spin valve signal much as the spin 
independent interface resistance just adds to the (larger) 
spin indep enden t resistance of the Cu strip of about 6 f2 
(see Sec. VIP eqs. |2^ and |9|). The spin signal can 
therefore only be significantly be reduced due to a pos- 
sible spin flip scattering mechanism at the interface, an 
effect yHHF nas recently been studied in CPP-GMR spin 
valves.P 8 lp 9 l The physical origin of this mechanism could 
be diverse, for instance: surface roughness creating local 
magnetic fields due to the formation of random domains, 
a dead magnetic layer of the first few nm of Co or Ni 
and the formation of anti-ferromagnetic oxides CoO and 
NiO at the surface during the time in between the Kauf- 
mann sputtering and the Cu deposition. However as we 
do not have a characterization of the interfacial struc- 
ture we cannot analyze what could be the most probable 
cause. We do note that the above mentioned mechanisms 
could also apply for the Py/Cu interface, however the 
Py/Cu/Py spin valve data show that their manifestation 
in these samples is apparently absent or less severe. 

Secondly, a change in the bulk properties of the Co 
and Ni could explain the small spin valve signals. All 
though in our opinion it is not likely that the bulk spin 
relaxation length would be subdue to a substantial short- 
ening, a reduction of the polarization qj? in our Co and Ni 
ferromagnets might occur. In CIP-GMR experimentsEll 
a strong decrease of more than an order of magnitude 
in the GMR signal was reported upon changing the base 
(H2O) pressure of the in the vacuum chamber from 10~ 8 
to 10~ 5 mbar, just before deposition the Co and Cu 
layers. In our deposition chamber the base pressure is 
only 10 mbar, whereas in the experiments e.g. on 
Co/Ag multilayerstil the base pressure fs-Qf the system 
is 10~ 8 mbar. However, theoretical workEZI does predict 
aco ~ 0.6 for Co with a conductivity close to our thin 
Co film a Co = 7.3 • 10 6 Qr x m~ x . We note howevec,that 
the Co layers in the Co/Cu multilayered nanowiresE 2 ! and 
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the Co/Ag multilayers^!! have a hep structure, whereas 
the calculations of Ref. ^ were done on fee cobalt. Un- 
fortunately we do not know the crystallinity and/or the 
crystal orientation of our Co films. 

X. CONCLUSIONS 

We have demonstrated spin injection and accumula- 
tion in metallic mesoscopic spin valves with transparent 
contacts. We have shown that in a conventional mea- 
surement geometry the magneto resistance effects of the 
injecting and detecting contacts can be much larger than 
the spin valve effect, making it impossible to observe the 
spin valve effect in a 'conventional' measurement geome- 
try. However, these contact effects can be used to moni- 
tor the magnetization reversal process of the spin inject- 
ing and detecting contacts. In a non-local measurement 
geometry we can completely isolate the spin valve effect, 
as was reported earlier in Ref. [l8| Using this geometry 
we find spin relaxation lengths in Cu of around 1 fim at 
T = 4.2 K and 350 nm at RT and spin relaxation lengths 



in Al of around 1.2 fim at T = 4.2 K and 600 nm at RT. 
The associated spin relaxation times in Al and Cu are in 
good agreement with theory and values from experiments 
previously reported in the literature. For the Py material 
we find spin relaxation lengths and current polarizations 
in agreement with GMR experiments. However for Co we 
obtain values of cif^f which are up to a factor 40 smaller 
than their GMR counterpart. For Ni electrodes we are 
unable to resolve a spin valve signal within the limits of 
our experimental accuracy, corresponding with olf Af at 
least a factor 10 lower than expected. Finally, we believe 
that the use of tunnel barriers should make it possible to 
increase the polarization of the injected current is non- 
magnetic metals, as we recently have shownHJ'EJ This 
should make it possible to increase the spin signals to 
about 1 Q in metals. 
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